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The recently discovered chimera state involves the coexistence of synchronized and desynchronized states 
for a group of identical oscillators. This fascinating chimera state has until now been found only in non-local or 
globally coupled oscillator systems. In this work, we for the first time show numerical evidence of the existence 
of spiral wave chimeras in reaction-diffusion systems where each element is locally coupled by diffusion. This 
spiral wave chimera rotates inwardly, i.e., coherent waves propagate toward the phase randomized core. A 
continuous transition from spiral waves with smooth core to spiral wave chimeras is found as we change the 
local dynamics of the system. Our findings on the spiral wave chimera in locally coupled oscillator systems 
largely improve our understanding of the chimera state and suggest that spiral chimera states may be found in 
natural systems which can be modeled by a set of oscillators indirectly coupled by a diffusive environment. 

PACS numbers: 05.45.Xt, 89.75.-k 


I. INTRODUCTION 

Collective behavior, which occurs commonly in physical, 
chemical and biological systems, has been a subject of contin¬ 
ued interest in nonlinear science over the last decades ffl. 
In neural and biological systems, a typical collective behav¬ 
ior observed is the coherent motion of oscillators. This phe¬ 
nomenon of synchronization has been widely regarded as hav¬ 
ing important implications to the function and performance of 
those systems jj, @0 . For instance, asynchronous contraction 
of the heart the may be triggered by electrical spiral waves, 
which eventually leads to heart dysfunction 0. 

Recently, much attention has been paid to a particu¬ 
lar hybrid state in which an ensemble of identical oscilla¬ 
tors with identical coupling spontaneously degenerate to one 
group with synchronization (coherent) and the other group 
with desynchronization (incoherent) IB8UT8n . This fascinating 
counterintuitive state was first discovered by Kuramoto and 
co-workers d], and named “chimera state” by Strogatz 19]. 
The experimental confirmation of the existence of chimera 
states has been reported independently almost at the same 
time in diverse systems such as coupled maps [ 10], chemi¬ 
cal oscillators (Til 12] and mechanical pendulums [0]. In 
two-dimensional systems, chimera states take the form of spi¬ 
ral waves GMi, called spiral wave chimeras with phase¬ 
locking oscillators in the spiral arm but a phase-randomized 
spiral core Hi. An analytical solution for a spiral wave 
chimera was further demonstrated by Strogatz et al. using 
a Kuramoto-type phase equation with nonlocal coupling lfl7ll . 
Recently, the spiral wave chimera state has been reported ex¬ 
perimentally in chemical oscillators Il2ll and numerically in 
complex and even chaotic oscillators fl8ll where nonlocal cou¬ 
pling is introduced. 

It was initially believed that the chimera state arises from 
so-called nonlocal coupling, i.e. each oscillator in the sys¬ 
tem will be affected instantaneously by a group of oscillators 
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within certain interacted range. Therefore, nonlocal coupling 
is intermediate between the cases of local and global coupling. 
However, some recent works show that the non-locality con¬ 
ditions for occurrence of the chimera state can be further re¬ 
laxed m. For example, Sethia et al. have demonstrated a 
generalized chimera state called “amplitude mediated chimera 
state” in a population of oscillators even in the case of global 
coupling rather than the nonlocal coupling Q. 

It is worth pointing out, on one hand, that in the seminal 
work of Kuramoto et al ., a key assumption to observe spi¬ 
ral wave chimeras is that the third component of reaction- 
diffusion (RD) changes so fast that it can be eliminated adia- 
batically [ 16]. In this way, the three-component RD system is 
essentially reduced to an effective two-component one with an 
extra term represented by a nonlocal coupling. On the other 
hand, to observe chimera state in experiments, the realization 
of nonlocal coupling strongly relies on a computer G3Q. 
Therefore, it remains an open question whether spiral wave 
chimeras exist in more natural systems. In such systems, the 
time scale of each component in the system may be compara¬ 
ble and the coupling will be mediated by a natural law such as 
diffusion. Further more, spiral waves can rotate outwardly or 
inwardly |2j3], however all the spiral wave chimeras reported 
previously are outward. There is very little information on the 
existence of inwardly rotating spiral wave chimera. 

In this work, we report the first numerical evidence of the 
existence of spiral wave chimeras in three-component RD sys¬ 
tem where each oscillatory element is locally coupled by dif¬ 
fusion, i.e., by nearest neighbor interaction. Differing from 
previously observed spiral wave chimeras in nonlocally cou¬ 
pled systems, our study shows remarkable novel properties of 
spiral wave chimeras: first the observed spiral wave chimeras 
rotate inwardly, i.e., a coherent wave propagates toward the 
phase-randomized center. Second, the non-diffusing compo¬ 
nents in the RD system exhibit chimera properties, while the 
diffusing components show a coherent spiral structure. A con¬ 
tinuous transition from spiral wave with a smooth core to spi¬ 
ral wave chimeras is also identified. Our findings in RD sys¬ 
tems show that the occurrence of the chimera states does not 
require nonlocal or global coupling and therefore provide key 
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hints to explore the chimera state in the natural world. 


II. REACTION-DIFFUSION MODEL AND METHODS 

Model description. Our starting point is a three-component 
RD system in two spatial dimensions EH0: 

d t X = ct>{aX -aX 3 - bY -cZ) + D x V 2 X, 
d t Y = (X-Y), (1) 

rd t Z = (t>€ 2 (X-Z) + D z X 2 Z. 

These equations describe the evolution of concentrations of 
chemical reactants X , Y and Z , where X is an activator and 
Y and Z are inhibitors. Dx and Dz denote the diffusion 
coefficients of chemical species of X and Z, respectively. The 
dimensionless parameter r represents the characteristic time 
scale of the system variable Z, and will be assigned a finite 
non-zero value. In principle, the parameter can be absorbed 
to a, b , c, a , e\ and 62 which represent other parameters, but 
we still write it explicitly here to keep the same form as in 
Ref. ED]. This three-component RD system actually is an 
extension of the two-component FitzHugh-Nagumo model by 
coupling to a third variable Z, which was proposed to study 
pattern formation in the Belousov-Zhabotinsky (BZ) reaction 
dispersed in water droplets of a water-in-oil aerosol OT (AOT) 
microemulsion system (BZ-AOT system) El and to model 
spot dynamics in gas discharges [2j|]. 

This study will differ from previous work on pattern forma¬ 
tion in RD systems in two aspects. The first modification is 
related to the ratio of the diffusion coefficients S = Dx/Dz . 
Traditionally 5 > 1, but in our study we take S <C 1 or 5 = 0. 
Therefore, we consider Eqs. 0 from a dynamical point as a 
model describing a large number of oscillators that communi¬ 
cate with each other via a diffusive environment. Such kind 
of models may be related to various systems such as bacte¬ 
ria, yeast cells and chemical oscillators. The second impor¬ 
tant thing is on the time scale of r. Our second modifica¬ 
tion is to omit the key assumption that the third variable Z 
changes fast, i.e., r 0 Ily,i22]. Hence, in our case the 
above three-component RD system can no longer be reduced 
a two-component system with traditional nonlocal coupling. 

Numerical methods. We numerically integrated Eqs. 0 
in time using the explicit forward Euler method with spatial 
step Ax = Ay = 0.2 and time step At = (Ax) 2 /(5 Dz), in 
a grid of N x = 512 by N y = 512 oscillators. We use a five- 
point stencil to evaluate the Laplacian term in (Q]). Note that 
the diffusion term in dTJ is effectively implemented as the in¬ 
teraction of each oscillator with its four nearest neighbours on 
a rectangular lattice. In this paper, we focus on the dynamics 
of system as we vary the model parameters a, 62, with fixed 
value of (/) = 0.62, b = 3.0, c = 3.5, e x = 1.0, a = 4.0/3.0, 
r = 1.0. In the case of S = 0, changing the diffusion coeffi¬ 
cient Dz is equivalently to rescale the size of the system, and 
thus we keep the same value Dz = 0.5 through the work. 

Order parameter. To analyze the state of the coupled os¬ 
cillators, we introduce the oscillator phase 9 by tan 0(f) = 
Y{r)/X(r). Since several observed states will be incoherent 
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FIG. 1: (color online) Spiral wave chimeras in a three-component 
RD system 0 with a = 1.50 and €2 = 0.15 . (a) denotes snapshot 
of the variable X ; (b) is the magnification the square region in (a); 

(c) represents the variable of Z; (d) is the snapshot of phase defined 
tan 0(f) = Y(r)/X(r) and (e) is a cross section along the line in 

(d) by defining 4> = 0/+i,j — 0j,j. It is noted that spiral core is 
incoherent for the variable X, while it is coherent for the variable Z. 
(f) is the spatiotemporal pattern of X along the middle horizontal line 
of (a) showing wave propagation towards the spiral chimera core. In 
panels (a-d), the wave pattern rotates clockwise. 


at the level of the simulation grid, we consider the domain as 
a discrete set of oscillators at positions x = I Ax, y = J Ay 
with phase Opj = arctan(l/ 5 j/X/ 5 j). 

To quantitatively study the size of the region with incoher¬ 
ent oscillations, we define the spatiotemporal order parameter 
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where Opj is the oscillator phase as defined above. The 
notation (/, J) denotes the set of nearest neighbors include 
itself and 1/(2 m + 1) is a normalization factor where m 
denotes the number of oscillators with the nearest coupling 
along a given spatial dimension. In the present case, we 
set m = 2, since we know the diffusion term in Eqs. CO) 
was computed using a 5-point discrete Laplacian. Finally, 
the time-averaged order parameter was computed as 07 5 j = 

limAT-^oo xt f£ +AT °i,j(f)W> w ^re AT is a time in¬ 
terval, e.g, AT = 20000 in our simulations. This quan¬ 
tity was used to estimate the radius of the incoherent state 
R = (R x + Ry)/ 2 where R x ( R y ) denotes maxim distance 
between the grid points along the x (y) line when 07, j < 0.9. 
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FIG. 2: (color online) Transition of chimera spiral waves as a func¬ 
tion of a. In the case of small a = 0.90, a coherent spiral wave is 
observed (a) and it is gradually changed to a spiral wave chimera (b) 
as we increase the system parameter a = 1.40. Further increasing a 
led to the formation of completely disorder state (c) where a = 1.70. 
We here also show the corresponding time averaged order parame¬ 
ter cr(refer to (d-f)) and the order parameter but along the center of 
medium (g-i). Except of a, other parameters are taken as in Fig. 1. 


III. RESULTS 

Existence of chimera state in a RD system with inhibitor dif¬ 
fusion only (5 = 0). Figure 1 demonstrates the existence of 
a spiral wave chimera state in the locally coupled RD system 
o. The snapshot of variable X and its magnification around 
the core are shown in Fig. 1(a) and (b), respectively at time 
t = 10 5 after the initiation of the spiral wave. Note that the 
X-variable around the spiral core is discontinuous while it is 
smooth far from the core, which is similar to previously re¬ 
ported spiral waves chimeras in nonlocal systems. This prop¬ 
erty is also true for the other non-diffusing variable Y. (Figure 
not shown). 

However, in contrast to X and Y, for the diffusing inhibitor 
Z, the whole spiral pattern is smooth even in the region closed 
to the spiral core as shown in Fig. 1(c). This property is differ¬ 
ent from the spiral wave chimeras found in nonlocally coupled 
oscillator systems where all of variables describing the oscil¬ 
lator demonstrate a similar chimera characteristic. 

The corresponding phase distribution of Opj and the de¬ 
pendency of <F/,j = 6i,j — Oi+i,j on the position along 
a horizontal cross section through the center of the medium 
(J = Ny/2 ) is illustrated in Fig. l(d-e). As one notes that 
close to spiral center, the quantity <F/ 5 j are nonzero, which 
means the phases of the oscillators around the spiral core 
are discontinuous or incoherent; while far way from it, the 
phases becomes continuous or coherent as indicated by the 


zero quantity of j. This behavior is the defining property 
of a chimera state for a group of identical oscillators: some os¬ 
cillate in a coherent way while a localized subgroup oscillate 
incoherently H[]. 

It is worth pointing out that the spiral wave chimera shown 
in Fig. 1 rotates inwardly, i.e. one sees the coherent waves 
propagate toward the phase-randomized core. This can also 
be seen from the spatiotemporal pattern for X variable along 
a horizontal cross section through the center of the medium 
(J = Ny/2) shown in Fig. 1(f). The inwardly rotating spi¬ 
ral waves, also called antispiral waves, have been reported 15 
years ago Ei. However, to our best knowledge, this is the 
first time to report such inwardly rotating wave chimeras in a 
RD system. 

Impact of local dynamics on the formation of spiral wave 
chimera. To get more insights to the formation and robust¬ 
ness of the spiral wave chimera, we investigate its behavior 
under a sweep of the parameter a. From 12111 . we recall that 
the system (Q} changes from a stationary state to an oscillatory 
state via a Hopf bifurcation when the value of the parameter 
a is increased. For a < 1, we observe spiral waves (inward) 
with a smooth core. From, a < 1.1 onwards, spiral chimera 
states are found, whose core radius increases with increasing 
a; compare panels (a) and (b) from Fig. 2. The incoherent re¬ 
gion of the chimera state monotonically grows with increasing 
a, until the incoherent state fully occupies the domain. This 
completely incoherent state is already reached in Fig. 2(c) 
for a = 1.7. We have found continuous transition from nor¬ 
mal spiral waves to spiral wave chimeras by controlling the 
parameter a. 

In Fig. 2(d-f), we show the time-averaged order parameter 
for different values of a. It is noted that in the incoherent 
chimera core, cfpj is smaller than one whereas in the outer 
region ap j « 1 due to the continuity of phase. Therefore, in 
the case of spiral waves, the region with apj ~ 1 is extremely 
small, see Fig. 2(d). On the contrary, for the turbulent state, 
(f/ 5 j < 1 everywhere. These statements can also be seen 
clearly from Fig. 2(g-i), which show the local order parameter 
apj along the line J = N y /2. 

The dependency of the chimera core radius on the param¬ 
eter a is shown in Fig. 3. When 0.90 < a < 1.0, the ra¬ 
dius R « Ax and thus in such a case, a classical spiral wave 
emerges. However, when a lies between 1.1 and 1.5, the ra¬ 
dius of incoherent spiral core is clearly finite, corresponding 
to the spiral wave chimera state. Further, we observe a fast 
transition from spiral waves chimera to completely turbulent 
state beyond a = 1.5. 

To give the global picture about how the local dynamics pa¬ 
rameter affects the behaviors of the spiral wave chimeras, we 
further identified the occurrence of the chimera state in a wide 
range of the a-e 2 parameter space. The results are summa¬ 
rized in Fig. 4. In this figure, full circles represent the stable 
spiral wave chimeras, the full triangles and squares denote the 
spiral waves and incoherent state. The cross means that os¬ 
cillations are not sustained. In the parameter space, the sep¬ 
aration line between spiral and spiral wave chimera is almost 
vertical, which means that a plays the key role in determining 
such state, while 62 has a more important role in determining 
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FIG. 3: (color online) Radius of the incoherent chimera core as a 
function of model parameter a, for 62 = 0.15. 
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FIG. 4: (color online) Phase diagram for the spiral wave chimera 
state in a — 62 parameter space showing occurrence of spiral 
waves(A), fully incoherent state (□), spiral wave chimera states 
(•) and a stable resting state (x). All the pattern were observed at 
t — 10 5 after the initiation of the spiral wave. 


whether spiral chimeras or turbulent states appear. 

Spiral wave chimeras with finite activator diffusion (S). 
Till now, we only consider the RD systems where only the 
diffusion of inhibitor, Z, is presented, i.e., 5 = 0. A natu¬ 
ral question is: how such kind of spiral wave chimeras would 
change when we increase S from zero to a finite value? Typi¬ 
cal results are presented in Fig. 5. For 62 = 0.45, spiral wave 
chimeras survive as we increase to S = 10 -3 , see Fig. 5(a) 
and magnification of the core region 5(d). The core region be¬ 
comes smooth but somehow spilt as we increase to S = 10 -2 , 
see Fig. 5(b) and 5(e). Further increasing S = 10 “ 1 finally 
leads to the appearance of the spiral wave with smooth core 
as illustrated in Fig. 5(c) and 5(f). These results are within 
our expectation as the diffusion coupling of the activator has 
a clear smoothing effect. 


FIG. 5: (color online) The impact of finite S on spiral wave chimera, 
(a) S = 10 -3 Spiral wave chimeras; (b) S = 10 “ 2 Spiral wave with 
split core; (c) 6 = 10 “ 1 Spiral wave with smooth core, (d-e) are the 
corresponding core region as indicated by the white square region in 
(a-c). e 2 = 0.45, a = 1.70. 


IV. DISCUSSION 

In this study, we have found states in locally coupled RD 
systems that are highly similar to the chimera states found in 
non-locally coupled oscillator systems id. We have identi¬ 
fied this spiral wave chimera state for a broad parameter range 
in our model that is based on a chemical system. Our re¬ 
sults are not limited to the particular RD system (Q]), when 
it is are linearly coupled to the a third variable that diffuses. 
Since no parameter tuning was involved in our numerical ex¬ 
periments, we believe that the robust nature of spiral wave 
chimeras found in locally coupled systems indicates that spiral 
wave chimeras may be a frequent self-organized spatiotempo- 
ral pattern in complex systems, which may arise from weak 
coupling of a set of oscillators by a diffusing variable. There¬ 
fore, we expect that natural chimera states will be found in 
oscillatory biological or chemical systems. 

Previous works on chimera states were focused on nonlo¬ 
cal or global coupling. In both cases, each element in the 
system was treated as a discrete oscillator. Whether chimera 
states exist in a spatially extended continuous media remains 
an open question. In a recent note d, one-dimensional 
chimera states in networks of pure local coupling were re¬ 
cently reported. However, this author did not consider spiral 
chimera states. 

In our work, we started from the spatial discretization of a 
Laplacian operator in Eq. (Q]) on a Cartesian grid, and found 
that the resulting wave became incoherent at the level of the 
grid scale. Furthermore, we have verified that using a finer 
discretization, the chimera state still persists. Thus, in locally 
coupled systems, the elusive chimera state breaks the contin¬ 
uum representation of the medium and forces one to consider 
the domain as a set of individual oscillators. Therefore, future 
analytical approaches will need to match the continuum dy¬ 
namics in the outer region with discrete non-linear dynamics 
in the spiral core region. 

It is noted that although we take the form of local coupling 
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between the Z variable, in the limit case of 6 = 0, such kind of 
local coupling of the inhibitor Z may give rise to the nonlocal 
effect due to the in absence of diffusion of activator X. How¬ 
ever, these nonlocal effects, which largely differs from the tra¬ 
ditional nonlocal coupling (where each oscillator in the sys¬ 
tem will be affected instantaneously by a group of oscillators 
within certain interacted range), are not only spatial but also 
temporal, because of the finite r. The systems with the local 
coupling that nevertheless demonstrate spatiotemporally non¬ 
local effects have been largely overlooked in the last decades. 
However, such systems may be very common in biological 
systems and thus deserve further investigation in the future. 

There may be various reasons for not observing chimeras 
in natural experimental settings before. From our study we 
note two necessary conditions: first, the local dynamics of 
each element needs to be oscillatory. Secondly, one requires 
close-to-vanishing spatial coupling of the observed variable. 

Finally, we note the essential property of the system (1) 
used in our study is that each oscillator (described by X and 
Y) is coupled indirectly via a nonuniform dynamical environ¬ 
ment which described by the variable Z. The systems with 
such property may denote a broad class of the systems such 
as chemical oscillators e.g., BZ-AOT [20] and BZ particles 
imersed in the solutions [3h engineered gene network [ 25th 
yeast cells [26]. From the results in this study, we expect that 
natural chimera states can be found in many oscillatory bio¬ 
logical or chemical systems. 


V. CONCLUSIONS 

In summary, we have shown that (inward) spiral wave 
chimeras do exist in spatially extended oscillatory media 
where only nearest-nearest interaction between the element is 
present. In such system, the nondiffusing components appear 
as a spiral wave chimera, while the diffusing variables show a 
coherent spiral wave structure. 

A continuous transition from coherent spiral waves to spiral 
waves chimeras is observed as we increase the model param¬ 
eter a that controls the Hopf bifurcation. A phase diagram 
for spiral wave chimeras is identified in the wide parameter 
space. We further discussed the smooth effects on the spi¬ 
ral core by the diffusion of the activator X. Our results on 
the spiral wave chimera in locally coupled oscillator systems 
largely improves our understanding of the chimera state and 
provides indications that the chimera state can be found in nat¬ 
ural systems such as biological and chemical systems where 
each oscillating element communicates via a diffusive dynam¬ 
ical environment. 
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